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ABSTRACT 

The "Mouse" (PWN G359.23— 0.82) is a spectacular bow shock pulsar wind nebula, powered by 
the radio pulsar J 1747— 2958. The pulsar and its nebula are presumed to have a high space velocity, 
but their proper motions have not been directly measured. Here we present 8.5 GHz interferometric 
observations of the Mouse nebula with the Very Large Array, spanning a time baseline of 12 yr. We 
measure eastward proper motion for PWN G359.23— 0.82 (and hence indirectly for PSR J1747— 2958) 
of 12.9 ± 1.8 mas yr^ 1 , which at an assumed distance of 5 kpc corresponds to a transverse space 
velocity of 306 ± 43 km s _1 . Considering pressure balance at the apex of the bow shock, we calculate 
an in situ hydrogen number density of approximately 1.0 ~to '% cm~ 3 for the interstellar medium through 

which the system is traveling. A lower age limit for PSR J1747— 2958 of 163i?,Q kyr ^ s calculated by 
considering its potential birth site. The large discrepancy with the pulsar's spin-down age of 25 kyr 
is possibly explained by surface dipole magnetic field growth on a timescale «15 kyr, suggesting 
possible future evolution of PSR J1747— 2958 to a different class of neutron star. We also argue that 
the adjacent supernova remnant G359. 1—0.5 is not physically associated with the Mouse system but 
is rather an unrelated object along the line of sight. 

Subject headings: ISM: individual (G359. 1—0.5, G359.23— 0.82) pulsars: individual 

(PSR J1747— 2958) — stars: neutron — supernova remnants 



1. INTRODUCTION 

The evolution of neutron stars and the potential rela- 
tionships between some of their observed classes remain 
outstanding problems in astrophysics. Proper motion 
studies of neutron stars can provide independent age es- 
timates with which to shed light on these questions. In 
particular, the well defin ed geometry of bow sho ck pul- 
sar wind nebulae fPWNe: IGaensler "fc Slanc 2006), where 
the relativistic wind from a high-velocity pulsar is con- 
fined by ram pressure, can be used as a probe to aid in 
the understanding of both neutron star evolution and the 
properties of the local medium through which these stars 
travel. 

The "Mouse" (PWN G359.23-0.82), a non-thermal 
radio nebula, was discovered as part of a ra- 
dio continuum survey of the Galactic center region 
(jYusef-Zadeh fc Ballvl Il987h . and was suggested to be 
owered by a young puls ar following X-ray detection 
Predehl fe Kulkarni| 119951 ). It is now recognized as a 



bow shock PWN moving supersonically throu gh the in- 
terstellar medium (ISM: IGaensler et al"1l2004D . Its ax- 
ially symmetric morphology, shown in Figure [TJ con- 
sists of a compact "head", a fainter "body" extend- 
ing for ~10", and a long "tail" that extends west- 
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ward behind the Mouse for ~40" and ~12' at X-ray 
and radio wavele ngths respectively (jGaensler et al.ll2004t 
iMori et aT1l2005l ). The cometary tail appears to indicate 
motion away from a nearby supernova remnant (SNR), 
G359. 1-0.5 (|Yusef-Zadeh fc Ballvlfl987l) . 

A radio pulsar, J1747— 2958, has been discovere d 
within the "head" of the Mouse (jCamilo et all [20021. 
PSR J1747-2958 has a spin period P = 98.8 ms and 
period derivative P = 6.1 x 10~ 14 , implying a spin-down 
luminosity E = 2.5 x 10 36 ergs s -1 , surface dipole mag- 
netic field strength B = 2.5 x 10 12 G, and characteristic 
age t c = P/2P = 25 kyr (iCamilo et al.l 120021 s ee also 
updated timing data from IGaensler et al.l I2004H . The 
distance to the pulsa r is >4 kpc from X-ray absorption 
(IGaensler et al.ll2004D . and <5.5 kpc from HI absorption 
(|Uchida et alj [l992). Here we assume that the system 
lies at a distance of d = 5d$ kpc, where dg = 1 ± 0.2 
(la). 

Given such a small characteristic age, it is natural to 
ask where PSR J1747— 2958 was born and to try and 
find an associated SNR. While it is possible that no 
shell -type SNR is visible, s uch as with the Crab pul- 
sar (iSankrit fc Hesterl 119971 ) and other young pulsars 
(jBraun et al.l 119891 ) . an association with the adjacent 
SNR G359. 1—0.5 appears plausible. This remnant was 
initially suggested to be an unrelated backgrou nd ob- 
ject near the Galactic center (lUchida et al.|[j~992l ). How- 
ever, it is now believed that the two may be located 
at ro ughly the same distance (jYusef-Zadeh fc Gaensler! 
120051 and references therein). By determining a proper 
motion for PSR J1747— 2958, this association can be sub- 
jected to further scrutiny (for example, see analysis of 
PSR B1757-24, PWN G5. 27-0.90 a nd SNR G5.4-1.2; 
Blaz ek et al.|[200l iZeiger et al.l [20081 ) . 
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Fig. 1.— VLA image of the Mouse (PWN G359.23-0.82) at 1.4 GHz with a resolution of 12'.'8x8'.'4 (reproduced from lGaensler et al.l 
120041) . The brightness scale is logarithmic, ranging between —2.0 and +87.6 mjy beam -1 as indicated by the scale bar to the right of the 
image. The eastern rim of SNR G359. 1-0.5 is faintly visible west of ~RA 17 h 46 m 25 s . 



As PSR J1747— 2958 is a very faint radio source, it 
is difficult to measure its proper motion interferometri- 
cally. It is also difficult to use pulsar timing to measure 
its proper motion due to timing noise and its location 
near the ecliptic plane (|Camilo et al. I [20021 ). To circum- 
vent these issues, in this paper we investigate dual-epoch 
high-resolution radio observations of the Mouse nebula, 
spanning 12 years from 1993 to 2005, with the intention 
of indirectly inferring the motion of PSR J1747-2958 
through the motion of its bow shock PWN. In § [2] wc 
present these observations. In §[3] we present our analysis 
and measurement of proper motion using derivative im- 
ages of PWN G359.23-0.82. In §g]we use our measure- 
ment to determine an in situ hydrogen number density 
for the local ISM, to resolve the question of a possible as- 
sociation with SNR G359. 1—0.5, and to investigate the 
age and possible future evolution of PSR J1747— 2958. 
We summarize our conclusions in § [5] 

2. OBSERVATIONS 

PWN G359.23-0.82 was observed with the Very Large 
Array (VLA) on 1993 February 2 (Program AF245) and 
again 7 on 2005 January 22 (Program AG671). Each of 
these observations were carried out in the hybrid BnA 
configuration at a frequency near 8.5 GHz. The 1993 
and 2005 epochs used on-source observation times of 
3.12 and 2.72 hours respectively. The 1993 observa- 
tion only measured RR and LL circular polarization 
products, while the 2005 observation measured the cross 
terms RL and LR as well. Both observations used the 
same pointing center, located at RA = 17 h 47 m 15!764, 
Dec = -29°58'1'/12 (J2000), as well as the same pri- 
mary flux calibrator, 3C286. Both were phase-referenced 
to the extragalactic source TXS 1741—312, located at 
RA = 17 h 44 m 23!59, Dec = -31°16'35"97 (J2000), which 
is separated by 1?4 from the pointing center. 

Data reduction was carried out in near identical 
fashion for both epoc hs using the MIRIAD package 
(jSault fc Killeenl [2006) . taking into consideration the 
slightly different correlator mode used in the 1993 data. 
This process involved editing, calibrating, and imaging 
the data using multi-frequency synthesis and square pix- 
els of size 50 x 50 milli-arcseconds. These images were 

7 An observation of PWN G359.23-0.82 was also carried out 
on 1999 October 8 (Program AG571). However, the target was 
observed mainly at low elevation and the point spread function and 
spatial frequency coverage were both poor as a result, thus ruling 
out the observation's amenability to astrometric comparison. 



then deconvolved using a maximum entropy algorithm 
and smoothed to a common resolution with a circular 
Gaussian of full width at half-maximum (FWHM) Of! 81. 
The resulting images are shown in the left column of Fig- 
ure!! 

The peak flux densities of the 1993 and 2005 im- 
ages are 3.24 and 3.25 mjy beam -1 , respectively; the 
noise in these two images are 51 and 35 //Jy beam -1 , 
respectively. The pulsar J 1747— 2958 is located at 
RA = 17 h 47 m 15?882, Dec = -29°58T/0 (J2000), within 
the region of intens e synchrotron emissi on seen in each 
image (see § 3.5 of iGaensler et al.ll2004h . Qualitatively 
comparing each epoch from the left column of Figure 
it appears that the head of PWN G359.23-0.82 has the 
same overall shape in both images, with a quasi-parabolic 
eastern face, approximate axial symmetry along a hori- 
zontal axis through the centre of the nebula (although 
the position of peak intensity seems to shift slightly in 
declination), and a small extension to the west. By eye 
the PWN seems to be moving from west to east over 
time, in agreement with expectation from the cometary 
morphology seen in Figure [TJ Beyond any minor mor- 
phological changes seen between the images in the left 
column of Figure (2) the Mouse nebula seems to have ex- 
panded slightly over time. 

3. ANALYSIS 

To quantify any motion between epochs, an algorithm 
was developed to evaluate the cross-correlation coeffi- 
cient over a range of pixel shifts between images, es- 
sentially by producing a map of these coefficients. This 
algorithm ma de use of the Karma visualization package 
(|Goochlll996f) to impose accurate non-integer pixel shifts. 

We applied our algorithm to the image pair from the 
left column of Figure [2] to determine an image offset mea- 
surement. To check that this offset measurement would 
be robust against any possible nebular morphological 
change between epochs, we also applied our algorithm 
to the same image pair when both images were equally 
clipped at various flux density upper- level cutoffs. We 
found that the offset measurement was strongly depen- 
dent on the choice of flux density cutoff. Clearly such 
variation in the measured shift between epochs was not 
desirable, as selection of a final solution would have re- 
quired an arbitrary assumption about the appropriate 
level of flux density clipping to use. There was also 
no strong indication that the region of peak flux den- 
sity in each image coupled with the exact location of 
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Fig. 2. — Left column: VLA observations of the Mouse at 8.5 GHz over two epochs separated by 12 yr. Each image covers a 14'.'5x8"5 
field at a resolution of 0"81, as indicated by the circle at the bottom right of each panel. The brightness scale is linear, ranging between 
—0.23 and +3.5 mjy beam -1 . Right column: Spatial x-derivative images for 1993 (top) and 2005 (bottom), covering the same regions as 
images in the left column. The images are shown in absolute value to increase visual contrast using a linear brightness scale spanning zero 
and the largest magnitude derivative from either image. The box between ~RA 17^47 m 15?8— 16?5 indicates the eastern region extracted 
for cross-correlation. All panels: Contours are shown over each image at 30%, 50%, 75% and 90% of the peak flux within their respective 
column group. The dashed vertical line in each panel has been arbitrarily placed at a right ascension near the brightness peak in the 
top-right panel in order to determine if motion can be seen by eye between epochs. 



PSR J1747— 1958. In order to isolate the motion of the 
pulsar from as much nebular morphological change as 
possible, we focused on a different method involving the 
cross-correlation of spatial derivatives of the images from 
the left column of Figure [H 

As there is not enough information to solve for an inde- 
pendent Dec-shift, we will only focus on an RA-shift, and 
will assume that any Dec changes are primarily due to 
morphological evolution of the Mouse nebula. To justify 
this assumption, we note simplistically that the cometary 
tail of the Mouse in Figure [T] is oriented at <5° from the 
RA-axis, and thus estimate that any Dec motion con- 
tributes less than 10% to the total proper motion of the 
nebular system. The small angle also justifies our deci- 
sion not to calculate derivatives along orthogonal axes 
rotated against the RA-Dec coordinate system. 

For each epoch, an image of the first spatial deriva- 
tive of intensity in the x (RA) direction was created by 
shifting the original image by 0.1 pixels along the RA- 
axis, subtracting the original from this shifted version, 
and then dividing by the value of the shift. These x- 
derivative images are shown in the right column of Fig- 
ured! where the brighter pixels represent regions of larger 
x-derivative from the corresponding left column images 
(note that these derivative images are shown in absolute 
value so as to increase their visual contrast; this opera- 
tion was not applied to the analyzed data). 

The x-derivative images have signal-to-noise ratios 
~13, since progression to higher derivatives degrades sen- 
sitivity. As seen in the right column of Figure [21 the 
x-derivative images of the Mouse are divided into two 
isolated regions: an eastern forward region and a west- 
ern rear region. Derivatives in the eastern region are 
greater in magnitude than those in the western region. 

As we will justify in § 01 we propose that the east- 
ern region of the x-derivative images tracks the forward 
termination shock of the Mouse nebula, which in turn 



acts as a proxy for an upper limit on the motion of 
PSR J1747— 2958. The eastern region provides a natural 
localized feature at each epoch with which to generate 
cross-correlation maps in order to track the motion of 
PSR J1747-2958. 

3.1. Calculation of Nebular Proper Motion 

To prepare the x-derivative images for cross- 
correlation, their eastern regions, extending between 
Right Ascensions (J2000) of 17 h 47 m 15!8 and 17 h 47 m 16?5 
(see Figure [U , were extracted and padded along each 
side with 50 pixels (2500 mas) of value zero. These 
cropped and padded x-derivative images for the 1993 
and 2005 epochs were then cross-correlated with each 
other over a range of non-integer pixel shifts between 
-2500 and +2500 mas in both RA and Dec. The resul- 
tant 2005—1993 cross-correlation map, which indicates 
the shift required to make the 1993 epoch colocate with 
the 2005 epoch, is shown in the left of Figure [U 

Note that the map in Figure [3] incorporates trial shifts 
large enough to probe regions where the cross-correlation 
falls to zero (corresponding to cross-correlation between 
signal and a source-free region, as opposed to only prob- 
ing trial shifts close to the maxima of each x-derivative 
map) . In this way, the contours presented in Figure [3] 
represent percentages of the peak cross-correlation value. 

To quantify the shift between epochs, a profile along a 
line parallel to the RA-axis was taken through the peak of 
the cross-correlation map, as shown in the right of Figure 
OH We assume that morphological changes are negligible 
in the eastern region of the x-derivative images; therefore, 
by taking a profile through the peak we tolerate small 
Dec shifts between the two epochs. 

The RA shift between the 2005 and 1993 epochs was 
determined by fitting a Gaussian to the central 660 mas 
of the cross-correlation profile. The resultant shift is 
154 mas with a statistical uncertainty of 7 mas, where 
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Fig. 3. — Left: Cross-correlation map for cropped and padded x-derivative images between 1993 and 2005. Shifts range from —2500 
to +2500 mas in both RA and Dec. The pixels are scaled from zero (white) to the peak value of the map (black), and contours are 
at 50%, 68.3%, 99.5%, and 99.7% of this peak value. Right: Profile along a line parallel to the RA-axis through the peak value of the 
cross-correlation map (solid curve). The caption quantifies the RA-shift for the fitted peak value (also indicated by the vertical dashed 
line), FWHM (also indicated by the dotted vertical lines), and signal-to- noise ratio. 



the latter is equal to the FWHM divided by twice the 
signal-to-noise ratio. This calculation reflects the angu- 
lar resolution and noise in the images, but to completely 
quantify the error on the shift between the two epochs, 
systematic errors also need to be incorporated. 

To estimate the positional error in the image plane, 
corresponding to phase error in the spatial frequency 
plane, the phases of the complex visibility data for each 
epoch were self-calibrated 8 . By dividing the standard 
error of the mean of the phase variation in the gain solu- 
tions by 180°, the fraction of a diffraction- limited beam 
by which positions may have been in error in the im- 
age plane were calculated. By multiplying this fraction 
with the diffraction limited beamwidth for the 1993 and 
2005 epochs, the two-dimensional relative positional un- 
certainty of these two reference frames was estimated to 
be 22 and 20 mas respectively. 

The systematic error in our measurement, which de- 
scribes the relative positional error between the two 
epochs, was then determined by reducing the self- 
calibrated positional uncertainties for each epoch by y/2 
(we are only looking at random positional uncertainties 
projected along the RA-axis), and added in quadrature. 
This error was found to be 21 mas, which totally domi- 
nates the much smaller statistical error of 7 mas. 

By calculating the total error as the quadrature sum of 
the statistical and systematic errors, the RA-shift of the 
PWN tip between the 1993 and 2005 epochs was found to 
be 154 ±22 mas. When divided by the 4372 days elapsed 
between these two epochs, the measured shift corre- 
sponds to a proper motion of /x = 12.9 ±1.8 mas yr _1 
in an eastward direction. We therefore detect motion 
at the 7cr level. Note that, at an assumed distance of 
^5 kpc along a line of sight toward the Galactic center, 
this motion is unlikel y to be contaminated sig nificantly 
by Galactic rotation (OUing_& Mcrrificld 1993). 

If we simplistically compare the eastward component 
of the proper motion with the angle of the Mouse's 
cometary tail, as described earlier in § [3l we obtain a 
crude estimate of ~1 mas yr _1 for the northerly com- 
ponent of the nebula's proper motion. As this value is 

8 Note that self-calibration was not used in the general reduction 
process because it would have caused a degradation in relative po- 
sitional information between the final images, as the phases would 
no longer be tied to a secondary calibrator of known position. 



well within the error for the eastward motion, which is 
dominated by systematic effects, we feel that our ear- 
lier assumption of pure eastward motion in the presence 
of relative positional uncertainty between the 1993 and 
2005 reference frames is justified. 

4. DISCUSSION 

Bow shock PWNe have a double-shock structure con- 
sisting of an outer bow shock where the ambient ISM is 
collisionally excited, an inner termination shock at which 
the pulsar's relativistic wind is decelerated, and a con- 
tact discontinuity between these two shocks which marks 
the boundary bet ween shocked ISM and shocked pul- 
sar wind material dGaensler k. S lane 200(f). The outer 
bow shock may emit in Ha, though for systems such 
as PWN G359.23-0.82 with high levels of extinction 
along their line of sight, the detection of such emission 
would not be expected. The inner termination shock, 
which encloses the pulsar's relativistic wind, may emit 
synchrotron radiation detectable at radio/X-ray wave- 
lengths. It is expected that any synchrotron emission 
beyond the termination shoc k would be sharply bo unded 
by the contact discontinuity (Gaensler et al. 2004]). 

As mentioned in § [3j we suggest that the eastern re- 
gions of the x-derivative images from Figure[5]provide the 
best opportunity to track motion of PSR J1747— 2958, 
relatively independent of any morphological changes oc- 
curring in PWN G359.23— 0.82. Physically, these regions 
of greatest spatial derivative (along the RA-axis) might 
correspond to the vicinity of the termination shock apex, 
or possibly the contact discontinuity between the two 
forward shocks, where motion of the pulsar is causing 
confinement of its wind and where rapid changes in flux 
might be expected to occur over relatively small angu- 
lar scales. This is consistent with hydrodynamic sim- 
ulations which predict that the apex of the bow shock 
will be located just outside a region of i ntense syn- 
chrotron emission in which the p ulsar lies (|Bucciantinil 
2002? Ivan der Swaluw et al J 12003 ). 

The assumption that the eastern region of each x- 
derivative image can be used as a proxy to track the 
motion of PSR J1747— 2958 is therefore plausible, but 
difficult to completely justify. To show that motion cal- 
culated in this way provides an upper limit to the true 
motion of PSR J1747— 2958 we recall the overall morpho- 
logical change described at the end of § [2j namely that 
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the Mouse nebula has expanded with time between the 
1993 to 2005 epochs. This expansion suggests that the 
ISM density may be dropping, causing the termination 
shock to move further away from pulsar, so that any mo- 
tion calculated using the nebula may in fact overestimate 
th e motion of the pul sar (a similar argument was used 
by iBlazek et al.l 120061 in placing an upper limit on the 
motion of the PWN associated with PSR B1757-24). 
Such changes in density are to be expected as the nebula 
moves through inter stellar space, where like the spectac- 
ular Guitar nebula (|Chatterjee fc Cordesl 120041 ) motion 
may reveal small-scale inhomogeneities in the density of 
the ISM. We therefore assume that our measurement of 
proper motion from § 13.11 corresponds to an upper limit 
on the true proper motion of PSR J1747-2958. 

4.1. Space Velocity and Environment of 

PSR J 1747- 2958 

Using our proper motion result from § 13. H and the argu- 
ments for interpreting this motion as an upper limit from 
§ H the projected eastward velocity of PSR J1747-2958 
is inferred to be V PSR}± < (306 ± 43) d 5 km s -1 . Given 
that no motion along the line of sight or in Dec could 
be measured, we will assume that our estimate of Vp Sai i_ 
approximates the 3-dimensional space velocity V PSB . 

In a bow shock PWN, the pulsar's relativistic 
wind will be confined and balanced by ram pres- 
sure. Using our proper motion upper limit (assum- 
ing Vp SR ~ U PSRi ), the pressure balance relationship 9 

Upsr. = ZQbn^^d^ 1 km s _1 from § 4.4 of lGaensler et al.l 
(2004), and Monte Carlo simulation, we find an in situ 
hydrogen number density n (I-O+q^) ^s" 4 cm ~ 3 at 

68% confidence, or n , 95% ~ (l-0icu) ^5 4 cm_3 at 95% 
confidence. Our calculated density no implies a local 
sound speed of ~5 km s _1 , corresponding to motion 
through the warm phase of the ISM. 

Our space velocity for PSR J1747— 2958 is compara- 
ble with other pulsars tha t have observed bow shocks 
(|Chatteriee fc Corded I2O 02:) , and is consistent with the 
overall proj ected velocity distr i bution of the young pulsar 
popu lation (|Hobbs et al . 2005; F aucher-Giguere fc Kaspil 
120061 ). We note that iGaensler et all (|2004l ) estimated a 
proper motion and space velocity of «25 mas yr _1 and 
«600 km s _1 , respectively, which are a factor of two 
larger than the values determined in this paper. How- 
ever, by halving their assumed sound speed of 10 km s , 
their estimates of motion correspondingly halve. 

We now use our proper motion and hydrogen number 
density results to resolve the question of association be- 
tween PSR J1747-2958 and SNR G359. 1-0.5, and to 
investigate the age and possible future evolution of this 
pulsar. 

4.2. Association with SNR G359.1-0.5? 

If PSR J1747-2958 and the adjacent SNR G359. 1-0.5 
are associated and have a common progenitor, then an 
age estimate for the system that is independent of both 
distance and inclination effects is simply the time taken 

9 This relationship assumes a uniform density p with typical 
cosmic abundances, expressed as p = 1.37nomj^, where mu is 
the mass of a hydrogen atom and no is the number density of the 
ambient ISM. 



for the pulsar to traverse the eastward angular separation 
between the explosion site inferred from the SNR mor- 
phology and its current location, at its inferred eastward 
proper motion < /1. Assuming pulsar birth at the cen- 
ter of the SNR, the eastward angula r separation b e tween 
the center of SNR G359. 1-0.5 from lUchidalTaT] (jl992l ) 
and the l ocation of PSR J 1747- 2958 from the timing so- 
lution bv IGaensler et all (|2004l ) is found to be ~ 23', 
which would imply a system age of 9 j \x > 110 kyr. 

Given such a large age, and the unremarkable interstel- 
lar hydrogen number density at the (currently assumed) 
nearby Mouse (from § I4.ip , it would be difficult to argue 
why SNR G359. 1—0.5 has not dissipated and faded from 
view. Instead, SNR G359. 1—0.5 appears to be a mid- 
dle aged remnan t ^18 kyr old which continues to emit 
thermal X-rays (jBamba et al.1 [2000) . We conclude, in- 
dependent of distance estimates to either the pulsar or 
the SNR, that PSR J1747-2958 is moving too slowly 
to be physically associated with the relatively young 
SNR G359. 1-0.5. 

4.3. Age Estimate for PSR Jl 747- 2958 

Given that an association between the Mouse and 
SNR G359. 1—0.5 is unlikely, as outlined in §021 we now 
estimate the age of PSR J1747— 2958 assuming that it is 
unrelated to this SNR. 

As seen in Figure [T] there is a cometary tail of emis- 
sion extending around 12' (~17ds pc) westward of the 
Mouse, containing shocked pulsar wind material flowing 
back from the termin ation shock about PSR J1747— 2958 
(|Gaensler et al.|[2004l ). We begin by simplistically assum- 
ing that this pulsar was born at the tail's western tip. 
By dividing the tail length by the proper motion /1, we 
estimate an age of t tail « 56 kyr. Note that this age is 
independent of any distance estimates to the Mouse or of 
inclination effects. However, given that the tail appears 
to simply fade away rather than terminate suddenly, it 
is possible that tail could be much longer, and thus that 
the system could be much older (by considering the up- 
per limit arguments from § |4j this syst em age may be 
even g reater still). As discussed in § 4.7 of lGaensler et al.l 
(2004) , it is unlikely that the Mouse or its entire tail could 
still be located inside an unseen associated SNR, given 
that the tail is smooth and uninterrupted. In addition, 
the lack of a rejuvenated SNR shell anywhere along the 
length of the tail (or indeed beyond it), such as that seen, 
for example, in th e interaction between PSR B1951+32 
and SNR CTB80 (|van der Swaluw et alJl2?M ). supports 
the conclusion that the Mouse's tail is located entirely 
in the ISM. Therefore, the rim of the Mouse's unseen 
associated SNR must be located at a minimum angular 
separation of ~12' west of the Mouse's current location, 
implying that f taU is a lower limit on the time elapsed 
since PSR J1747— 2958 was in the vicinity of this rim. 

To estimate the total age of PSR J1747-2958 we 
thus need to incorporate the time taken for this pul- 
sar to escape its associated SNR and reach its current 
location, taking into account the continued expansion 
of the SNR following pulsar escape (which will sweep 
up, and therefore shorten, part of the Mouse's tail ini- 
tially located in t he ISM). Using Monte Carl o simula- 
tion and following Ivan der Swaluw et all (|2003l ). we find 
that the time taken for a pulsar, traveling with veloc- 
ity Vpspl± < (306 ±43)^5 km s _1 through a typical in- 
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terstellar environment with constant hydrogen number 
density no from § 14. 11 to escape from its SNR (while 
in the pressure-drive snowplow phase) of typical explo- 
sion energy ^10 51 ergs, and leave behind a 12' tail in 
the ISM which remains ahead of the expanding SNR, 
is 163t^ kyr at 68% confidence, or 163^° kyr at 95% 
confidence. Note that the errors quoted for this total 
time incorporate the error in [i and are only weakly de- 
pendent on the uncertainty in distance d$ (for compar- 
ison, when the distance to PSR J1747— 2958 is fixed at 
5 kpc, the 68% and 95% confidence intervals are reduced 
to 164;t^ kyr and 164+^ kyr, respectively). 

Assuming that PSR J1747— 2958 was created in such a 
supernova explosion, and noting that the pulsars's travel 
time in the ISM t tM is a lower limit (even without tak- 
ing into account the upper limit associated with [i, as 
described earlier in this section), we can thus establish a 
lower limit on the age of the pulsar of i total > 163^20 kyr 
(68% confidence). This lower limit is greater than 6 times 
the characteristic age r c = 25.5 kyr of PSR J1747-2958, 
which was derived from its measured pulse P and P 
(jCamilo et al.l f2002). suggesting that, within the context 
of the characteristic approximation, the spindown prop- 
erties of this pulsar deviate significantly from magnetic 
dipole braking (see § 14.4)) . 

Our result is similar to the age discrepancy pre- 
vious ly claimed for PSR B1757-24 (jGaensler fc Frail 
2000); however, ambiguity regarding association with 
SNR G5.4— 1.2 presents difficulties with this claim 
(jThorsett et al.ll2002t IZeiger et al.ll2008D . In comparison, 
given the relatively simple assumptions made in this pa- 
per, PSR J1747— 2958 arguably provides the most robust 
evidence to date that some pulsars may be much older 
than their characteristic age. We now discuss the poten- 
tial implications of this age discrepancy with regard to 
the future evolution of PSR J1747-2958. 

4.4. Possible Future Evolution of PSR J 1747- 2958 

Pulsars are assumed to slow down in their rotation ac- 
cording to the spin-down relationship ft = —K£l n , where 
fi is the pulsar's angular rotation frequency, fi is the an- 
gular frequency derivative, n is the braking index, and 
K is a positive constant that depends on the pulsar's 
moment of inertia and magnetic moment. By taking 
temporal derivatives of this spin-down relationship, the 
rate at which the characteristic age r c chang es with tim e 
can be expressed as dr c /dt=(n — l)/2 (e.g., Eyni[2004). 
Evaluating dr c /dt as the ratio between r c = 25.5 kyr 
and the lower age limit i total , we estimate a braking in- 
dex of n < 1.3 for PSR J1747-2958 (incorporating the 
error limits from t total does not significantly affect this 
value). Given that magnetic dipole braking corresponds 
to n = 3 (the value assumed when calculating r c ), the 
smaller braking index calculated here indicates that ei- 
ther some form of non-standard pulsar braking is taking 
place, or that standard magnetic dipole braking is oc- 
curring in the presence of evolutio n of the magnetic field 
or of the moment o f inertia (e.g., iBlandford fc Roman! 
[19881: ICamildll99fih . If we adopt a constant moment 
of inertia and assume standard electromagnetic dipole 
braking, then by performing a similar derivation from 
the spin-down relationship for the surface magnetic field, 
the magnetic field growth timescale can be expressed as 



B/ (dB/dt)=r c / (3 - n) (e.g., Evnel |200l . Evaluating 
this with our braking index, the magnetic field growth 
timescale 10 is estimated to be «15 kyr. 

It is interesting to note that the braking index inferred 
here for PSR J1747— 2958 is comparable to the value ob- 
tained from an estimate of P made for the Vela pul- 
sar B0833— 45, w hich was found to have n — 1.4 ± 0.2 
(|Lvne et al.lll996D . To investigate the possible future evo- 
lution of PSR J1747— 2958, we plot its implied trajectory 
(along with that of the Vela pulsar) across the P—P dia- 
gram, as shown in Figured] Note that the magnitude of 
each plotted vector indicates motion over a timescale of 
50 kyr assuming that P is constant, which is not true for 
constant n; however, the trend is apparent. The longer 
vector for the Vela pulsar simply indicates that it is brak- 
ing more rapidly than PSR J1747-2958. 




A B0B33-45 (Vela) 

H J0537-6910 



10 3 0.01 0.1 1 10 

Period (s) 

Fig. 4. — The P — P diagram for rotating neutron stars, where 
points indicate the period and period derivative for over 1600 rotat- 
ing stars where such measurements were availa ble (data obtained 
from the ATNF Pulsar Catalogue, version 1.34; Manchest er et al.l 
120051) . The periods and period derivatives of three neutron stars, 
B0833-45 (the Vela pulsar), J1747-2958 (the Mouse pulsar), and 
J0537— 6910 (a 16-ms pulsar in the Large Magellanic Cl oud) are 
indicated with symbols (see discussion at the end of § 14.41 1. as 
well as those of the known magnetars. Dashed lines indicate con- 
tours of constant t c , while solid lines indicate contours of constant 
surface magnetic field B. Points toward the bottom left are mil- 
lisecond (recycled) pulsars, those in the central concentration are 
middle— old age pulsars, and those to the top left are young pul- 
sars. The middle— old age pulsars are presumed to move down 
to the right along lines of constant magnetic field, corresponding 
to magnetic dipole braking, while some young pulsars such as the 
Vela pulsar have been found to be moving away from the main 
body of pulsars, up and to the right toward the re gion of the mag- 
netars. The plotted vectors for PSR B0833-45 IILvne! I2004TI and 
for PSR J1747— 2958 (this paper) represent inferred trajectories 
across the diagram over 50 kyr, assuming a constant braking index 
over this time period; the magnitude of each vector is proportional 
to the magnetic field growth timescale of the respective neutron 
star. 

10 As noted by ILvnd (2003) the magnetic field growth is not 
linear over time, as can be appreciated by taking an example of a 
pulsar with braking index n=l. 
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The plotted vectors for the Vela and Mouse pulsars 
both seem to point in the direction of the magnetars 
(high-energy neutron st ars for which i?>10 14 G; for a re- 
view of magnetars, see iWoods fc Thompson! 120061 ) . By 
extrapolating the trajectories in Figure 21 and assum- 
ing negligible magnetic field decay over time, it can be 
suggested that young, ener getic, rapidly spinning pulsars 
such as PSR J0537-6910 ([Marshall et aLlH99l . whose 
location in the P — P plane is shown with a star symbol, 
may evolve into objects like the Vel a or Mouse pulsars, 
which, as proposed by ILvnd ([2004) , may in turn con- 
tinue to undergo magnetic field growth until arriving in 
the parameter space of the magnetars. 

5. CONCLUSIONS 

We have investigated two epochs of intcrferometric 
data from the VLA spanning 12 years to indirectly infer a 
proper motion for the radio pulsar J1747— 2958 through 
observation of its bow shock PWN G359.23-0.82. 
Derivative images were used to highlight regions of rapid 
spatial variation in flux density within the original im- 
ages, corresponding to the vicinity of the forward termi- 
nation shock, thereby acting as a proxy for the motion 
of the pulsar. 

We measure an eastward proper motion for 
PWN G359. 23-0.82 of (J, = 12.9 ± 1.8 mas yr" 1 , 
and interpret this value as an upper limit on the motion 
of PSR J1747-2958. At this angular velocity, we 
argue that PSR J1747— 2958 is moving too slowly to be 
physically associated with the relatively young adjacent 
SNR G359. 1—0.5, independent of distance estimates to 
either object or of inclination effects. 

At a distance d — kpc, the proper mo- 

tion corresponds to a projected velocity of 
Vpsr,j_ < (306 ± 43) d$ km s — 1 , which is consistent 
with the projected velocity distribution for young 
pulsars. Combining the time taken for PSR J1747— 2958 
to traverse its smooth ~12' radio tail with the time to 
escape a typical SNR, we calculate a lower age limit for 
PSR J1747-2958 of t total > 163j|g kyr (68% confidence). 

The lower age limit t totM exceeds the characteristic age 
of PSR J1747— 2958 by more than a factor of 6, arguably 
providing the most robust evidence to date that some 
pulsars may be much older than their characteristic age. 
This age discrepancy for PSR J1747— 2958 suggests that 
the pulsar's spin rate is slowing with an estimated brak- 
ing index n < 1.3 and that its magnetic field is grow- 
ing on a timescale as 15 kyr. Such potential for mag- 



netic field growth in PSR J1747— 2958, in combination 
with other neutron stars that transcend their archety- 
pal categories such as PSR J1718— 3718, a radio pulsar 
with a magnetar-strength magnet ic field that does not 
exhib it magnetar-like emission (jKaspi fc M cLaughlin 
120051) . PSR J1846-0258, a rotation-po wered pulsar 
that exhibits magnetar-lik e behaviour (|Gavriil et al.l 
120081: lArchibald et al.ll2008f) . and magneta rs such as IE 
1547 .0—5408 that exhibit radio emission (jCamilo et al.l 
120081 and references therein), supports the notion that 
there may be evolutionary links between the rotation- 
powered and magnetar classes of neutron stars. How- 
ever, such a conclusion may be difficult to reconcile 
with evidence suggesting that magnetars are derived 
from more massive p r ogenitors than nor mal pulsars (e.g., 
IGaensler et all 2005: M uno et alJl2006h . If the massive 
progenitor hypothesis is correct, then this raises further 
questioning of whether, like the magnetars, there is any- 
thing special about the progenitor properties of neutron 
stars such as PSR J1747— 2958, or whether all rotation- 
powered pulsars exhibit similar magnetic field growth or 
even magnetar-like phases in their lifetimes. 

To constrain the motion of PSR J1747-2958 further, 
future observational epochs are desirable. It may be pos- 
sible to better constrain the motion and distance to this 
pulsar by interferometric astrometry with the next gen- 
eratio n of sensitive radio telescopes (e.g., ICordes et al.l 
2004!). High time resolution X-ray observations may also 
be useful to detect any magnetar-like behaviour from this 
rotation-powered radio pulsar. In general, more neutron 
star discoveries, as well as measured or inferred braking 
indices, may allow for a better understanding of possible 
neutron star evolution. 
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